Mon. Not. R. Astron. Soc. 000, 1-?? (0000) 



Printed 2 February 2008 



(MN I^TeX style file v2.2) 



A different approach for the estimation of Galactic model 
parameters 



S. Karaali, 1 * S. Bilir 1 , and E. Hamzaoglu 2 

1 Istanbul University Science Faculty, Department of Astronomy and Space Sciences, 34-119, University -Istanbul, Turkey 
. 2 Faculty of Engineering and Design, Istanbul Commerce University, Ragip Gumu§pala Caddesi No: 84, 34378 Eminonii- Istanbul, Turkey 

O • 

o : 

CM ■ 

1 2 February 2008 

ABSTRACT 



> 

<N 



We estimated the Galactic model parameters by means of a new approach based on 
the comparison of the observed space density functions per absolute magnitude interval 
with a unique density law for each population individually, and via the procedure in 
situ for the field SA 114 (a = 22 /l 40 m 00 s ,£ = OO^O'OO"; I = 68°. 15, b = -48°.38; 
4.239 square- degree; epoch 2000). The separation of stars into different populations has 
1^1. ■ been carried out by their spatial distribution. The new approach reveals that model 

1 parameters are absolute magnitude dependent. The scale height for thin disk decreases 

monotonously from absolutely bright (M(g )=5) to absolutely faint (M(g ) = 13) 
stars in a range 265-495 pc, but there is a discontunity at the absolute magnitude 
M(g ) = 10 where the squared secans hiperbolicus density law replaces the exponential 
one. The range of the scale-height for thick disk, dominant in the absolute magnitude 
interval 5 < M(g ) < 9, is less: 805-970 pc. The local space density for thick disk 
£^ ■ relative to thin disk decreases from 9.5% to 5.2% when one goes from the absolutely 

\ bright to faint magnitudes. Halo is dominant in three absolute magnitude intervals, 

i.e. 5 < M(g') < 6, 6 < M(g') < 7, and 7 < M(g') < 8 and the axial ratio for this 
component is almost the same for these intervals where c/a ~ 0.7. The same holds for 
the local space density relative to the thin disk with range (0.02-0.15)%. The model 
parameters estimated by comparison of the observed space density functions combined 
for three populations per absolute magnitude interval with the combined density laws 
agree with the cited values in the literature. Also each parameter is equal to at least one 
of the corresponding parameters estimated for different absolute magnitude intervals 
by the new approach. We argue that the most appropriate Galactic model parameters 
are those, that are magnitude dependent. 

Key words: Keywords: Technique: photometric-survey - Galaxy: stellar content - 
Method: data analysis 



1 INTRODUCTION 

For some years, there has been a conflict among the re- 
searchers about the history of our Galaxy. Yet, there is a 
large improvement about this topic since the poineering 
work of Eggen, Lynden-Bell, & Sandage (1962, hereafter 
ELS) who argued that the Galaxy collapsed in a free-fall 
time (~ 2 x 10 8 yr). Now, we know that the Galaxy col- 
lapsed in times of many Gyr (e.g. Yoshii & Saio 1979; Norris, 
Bessel, & Pickles 1985; Norris 1986; Sandage & Fouts 1987; 
Carney, Latham, & Laird 1990; Norris & Ryan 1991; Beers 
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& Sommer-Larsen 1995) and at least some of its compo- 
nents are formed from merger or accretion of numerous frag- 
ments, such as dwarf-type galaxies (cf. Searle & Zinn 1978; 
Freeman & Bland-Hawthorn 2002, and references therein). 
Also, the number of population components of the Galaxy 
increased from two to three, complicating interpretations of 
any dataset. The new component, 'thick disk', introduced by 
Gilmore & Reid (1983) in order to explain the observation 
that star counts towards the South Galactic Pole were not 
in agreement with a single disk, 'thin disk', component but 
rather could be nicely represented by two such components. 
The new component is discussed by Gilmore & Wyse (1985) 
and Wyse & Gilmore (1986). 

The researchers use different methods to determine the 
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parameters for three population components and try to in- 
terpret them in relation to the formation and evolution of 
the Galaxy. Among the parameters, the local density and 
the scale-height of thick disk are the ones whose numeri- 
cal values improved relative to the original ones claimed by 
Gilmore & Reid (1983). In fact, the researchers indicate ten- 
dency to increase the original local density of thick disk from 
2% to 10% relative to the total local density and to decrease 
its scale- height from the original value 1.45 kpc down to 0.65 
kpc (Chen et al. 2001). In some studies, the range of the pa- 
rameters is large especially for the thick disk. For example, 
Chen et al. (2001) and Siegel et al. (2002) give 6.5% - 13% 
and 6% - 10% for the relative local density for thick disk, 
respectively. Now a question arises: what is the reason of 
these large ranges in two recent works (and in some other 
works) where one expects the most improved numerical val- 
ues? This is the main topic of our paper. We argue that, 
though considerable improvements were achieved, we still 
couldn't chose the most appropriate procedure for the esti- 
mation of Galactic model parameters. In the present study, 
we show that, parameters of Galactic model are functions 
of absolute magnitude and if this is not taken into consid- 
eration, different parameters with large ranges are obtained 
and this can not be unavoided. 

In sections 2 and 3 different methods and density law 
forms are discussed. The procedure used in this study is 
given in section 4. Section 5 provides de-reddened apparent 
magnitude, absolute magnitude, distance and density func- 
tion determinations. Model parameter estimation by differ- 
ent procedures is given in section 6 and finally section 7 
provides discussion. 



2 THE METHODS 

The studies related to the Galactic structure are usually 
carried out by star counts. Direct comparison between the 
theoretical and observed space densities is also used as an- 
other method. The first method is based on the fundamental 
equation of stellar statistics (von Seelinger 1898) which may 
be written as follows: 

A(mv, Sb-v) — ^^Ai(m v , Sb-v) 

= MM,S)D t (r)r 2 dr (1) 

where A is the differential number of counts at any partic- 
ular magnitude and colour, Ai is the contribution to those 
counts from population z, Q is the solid angle observed, <J>i 
is the luminosity function of population i, and Di is the 
density distribution of population i as a function of abso- 
lute magnitude M and spectral type 5, and r is distance 
along the line of sight. Here, the number of counts is the 
sum over the stellar populations of the convolution of the 
luminosity and density distribution functions. It is stated 
by many authors (cf. Siegel et al. 2002), the noninvertabil- 
ity and the vagaries of solving the nonunique convolution 
by trial and error, limit the star counts and will be a weak 
tool for exploring the Galaxy. The large number of Galactic 
structure models derived from star count studies confirms 
the nonuniquness problem (Table 1). The most conspicuous 



point in Table 1 is the large range of thick disk parame- 
ters, indicating a less certain density law for this compo- 
nent. Whereas the thin disk parameters occupy a narrow 
range of values. For halo, the results from star count sur- 
veys cover almost the entire range of parameter space from 
flattened de Vaucouleurs spheroid (Wyse & Gilmore 1989, 
Larsen 1996) to perfectly spherical power-law distributions 
(Ng et al. 1997). 

There is not enough study in the literature, carried out 
by comparison of theoretical and observational space den- 
sities. The works of Basle group (del Rio & Fenkart 1987, 
and Fenkart & Karaali 1987), and recently, Phleps et al. 
(2000), Siegel et al. (2002), Karaali et al. (2003), and Du et 
al. (2003) can be mentioned as examples. Photometric par- 
allaxes provide direct evaluation of spatial densities. Hence, 
the observations can be translated into discrete density mea- 
surements at various points in the Galaxy, instead of trying 
to fit the structure of the Galaxy into the observed param- 
eter space of colours and magnitudes. 

Almost the same results are seen in several studies. 
This is not a surprise, since such studies have explored sim- 
ilar data sets with similar limitations and additionally they 
probe the same direction in the sky such as the Galactic 
poles. Most studies are based on investigation of one or a few 
fields in different directions. The deep fields are small with 
corresponding poor statistical weight and the large fields are 
limited with shallower depth which may not be able to probe 
the Galaxy at large distances. The works of Reid & Majew- 
ski (1993) and Reid et al. (1996) can be given as examples for 
the first category, while the one of Gilmore & Reid (1983) 
for the second category. Deep surveys based on the mul- 
tidirectional Hubble Space Telescope (HST) (Zheng et al. 
2001), has another limitation, i.e. star-galaxy separation be- 
comes difficult at faint magnitudes. There are few programs 
which survey the Galaxy in multiple directions such as the 
Basle Halo Program (cf. Buser, Rong, & Karaali 1999), the 
Besancon program (cf. Robin et al. 1996, 2003), the APS- 
POSS program (Larsen 1996), and recently the SDSS (Chen 
et al. 2001). 

Star count studies in a single direction can lead to de- 
generate solutions and surveys, limited with small areas at 
the Galactic poles are insensitive to radial terms in the pop- 
ulation distributions (Reid & Majewski 1993, Robin et al. 
1996, and Siegel et al. 2002). 



3 THE DENSITY LAW FORMS 

Disk structures are usually parameterised in cylindrical co- 
ordinates by radial and vertical exponentials, 

Di(x, z) = me- z/Hi e- (x - Ro)/hi (2) 

where z is the distance from Galactic plane, x is the planar 
distance from the Galactic center, Ro is the solar distance to 
the Galactic center (8.6 kpc), Hi and hi are the scale height 
and scale length respectively, and m is the normalized local 
density. The suffix i takes the values 1 and 2, as long as thin 
and thick disks are considered. A similar form uses the sech 2 
(or sech) function to parameterize the vertical distribution 
for thin disk, 
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Table 1. Previous Galactic models. Symbols: TN: thin disk, TK: thick disk, S: spheroid (halo), R e \ effective radius, c/a: axes ratio. The 
figures in the parentheses for Siegel et al. (2002) are the corrected values for binarism. (*) power-law index replacing R e . 



H (TN) 


h(TN) 


n (TK) 


H(TK) 


h(TK) 


n (S) 


Re(S) 


c/a 


Reference 


(pc) 


(kpc) 




(kpc) 


(kpc) 




(kpc) 






310-325 




0.0125-0.025 


1.92-2.39 










Yoshii 1982 


300 




0.02 


1.45 




0.0020 


3.0 


0.85 


Gilmore & Reid 1983 


325 




0.02 


1.3 




0.0020 


3.0 


0.85 


Gilmore 1984 


280 




0.0028 


1.9 




0.0012 






Tritton & Morton 1984 


200-475 




0.016 


1.18-2.21 




0.0016 




0.80 


Robin & Creze 1986 


300 




0.02 


1.0 




0.0010 




0.85 


del Rio & Fenkart 1987 


285 




0.015 


1.3-1.5 




0.0020 


2.36 


Flat 


Fenkart et al. 1987 


325 




0.0224 


0.95 




0.0010 


2.9 


0.90 


Yoshii et al. 1987 


249 




0.041 


1.0 




0.0020 


3.0 


0.85 


Kuijken & Gilmore 1989 


350 


3.8 


0.019 


0.9 


3.8 


0.0011 


2.7 


0.84 


Yamagata & Yoshii 1992 


290 






0.86 






4.0 




VonHippel & Bothun 1992 


325 




0.0225 


1.5 




0.0015 


3.5 


0.80 


Reid & Majewski 1993 


325 


3.2 


0.019 


0.98 


4.3 


0.0024 


3.3 


0.48 


Larsen 1996 


250-270 


2.5 


0.056 


0.76 


2.8 


0.0015 


2.44-2.75* 


0.60-0.85 


Robin et al. 1996, 2000 


290 


4.0 


0.059 


0.91 


3.0 


0.0005 


2.69 


0.84 


Buser et al. 1998, 1999 


240 


2.5 


0.061 


0.79 


2.8 






0.60-0.85 


Ojha et al. 1999 


330 


2.25 


0.065-0.13 


0.58-0.75 


3.5 


0.0013 




0.55 


Chen et al. 2001 


280(350) 


2-2.5 


0.06-0.10 


0.7-1.0 (0.9-1.2) 


3-4 


0.0015 




0.50-0.70 


Siegel et al. 2002 



Di(x,z) =n i sech 2 (z/H i )e- (x - Ro)/hi . (3) 

As the secans hyperbolicus is the sum of two exponentials, 
H i is not really a scale height, but has to be compared to 
Hi by multiplying it with arcsech(l/e) ~ 1.65745 : Hi = 
1.65745#-. 

The density law for spheroid component is parame- 
terised in different forms. The most common is the de Vau- 
couleurs (1948) spheroid used to describe the surface bright- 
ness profile of elliptical galaxies. This law has been depro- 
jected into three dimensions by Young (1976) as follows: 

Ds(R) = n s exp[-7.669(R/R e ) 1/4 ]/(R/R e ) - 875 (4) 

where R is the (uncorrected) Galactocentric distance in 
spherical coordinates, R e is the effective radius, and n s is 
the normalised local density. R has to be corrected for the 
axial ratio k = c/a, 

R=[x 2 + (z/k) 2 ] 1/2 (5) 
where, 

x = [R 2 + (z/tanb) 2 - 2R Q (z/ tan b) cos if (6) 

b and I being the Galactic latitude and longitude for the 
field under investigation. The form used by the Basle group 
is independent of effective radius but the distance of the Sun 
to the Galactic center: 

D B (R) = ns exp[10.093(l - R/R ) 1/4 ]/(R/R ) - 875 (7) 

and alternative formulation is the power- law, 

D s (R) = n s /(a n +R n ) (8) 

where a is the core radius. 

Equations (2) and (3) replaced by (9) and (10) respec- 
tively, as long as vertical direction is considered, where: 

Di(z)=me- X/Hi (9) 



D z (z) = msech 2 {-z/H'i) (10) 



4 THE PROCEDURE USED IN THIS WORK 

In this work, we compared the observed and the theoretical 
space densities per absolute magnitude interval in the verti- 
cal direction of the Galaxy for a large range of absolute mag- 
nitude interval, 4 < M(g ) < 13, down to the limiting mag- 
nitude g = 22. The procedure is similar to that of Phleps et 
al. (2000), however the approach is different. First, we sep- 
arated the stars into different populations by a slight modi- 
fication of the method given by Karaali (1994), i.e. we used 
the spatial distribution of stars as a function of both absolute 
magnitude and apparent magnitude, whereas Karaali pref- 
ered a unique distribution for stars of all apparent magni- 
tudes. Second, we derived model parameters for each popu- 
lation individually for each absolute magnitude interval and 
we observed their difference. Third, and finally, the model 
parameters are estimated by comparison of the observed ver- 
tical space densities with the combined density laws (equa- 
tions 7, 9, and 10) for stars of all populations. In the last 
process we obtained two sets of parameters, one for the ab- 
solute magnitude interval 5 < M(g ) < 10, and another one 
for 5 < M(g ) < 13. We notice that the behavior of stars 
with 10<M(g)<13is different. We must keep in mind 
that many researchers related to estimation of the model 
parameters are restricted with M(V) < 8 (cf. Robin et al. 
2003). Different behavior of stars with 10 < M(g ) < 13 may 
result in different values and large ranges may be expected 
for the model parameters based on star counts. 

5 THE DATA AND REDUCTIONS 
5.1 Observations 

Field SA 114 (a = 22^40 m 00 s ,£ = 00°00 / 00 // ; I = 68°. 15, 
b = -48°. 38; 4.239 square-degree; epoch 2000) was mea- 
sured by the Isaac Newton Telescope (INT) Wide Field 
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Figure 1. E(B — V) colour-excess contours for the field SA 114 
as a function of Galactic latitude and longitude. 
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Figure 2. Apparent magnitude histogram for all sources. Bright 
stars are effected by saturation. 



Camera (WFC) mounted at the prime focus (f/3) of the 
2.5m INT on La Palma, Canary Islands, during seven ob- 
serving runs, i.e. 1998 September 3, 1999 July 17-22, 1999 
August 18, 1999 October 9-10, 2000 June 25-30, 2000 Oc- 
tober 19, and 2000 November 15-22. The WFC consists of 
4 EVV42 CCDs, each containing 2k x 4k pixels. They are 
fitted in a L- shaped pattern which makes the camera 6k x 
6k, minus 2k x 2k corner. The WFC has 13.5/i pixels cor- 
responding to 0.33 arcsec pixel -1 at INT prime focus and 
each covers an area of 22.8 x 11.4 arcmin on the sky. A 
total of 0.29 square-degree is covered by the combined of 
four CCDs. With a typical seeing of 1.0-1.3 arcsec on the 
INT, point objects are well sampled, which allows accurate 
photometry. 

Observations were taken in five bands (u RGO , g , r , i , 
z rgo) with a single exposure of 600s to nominal 5cr limiting 
magnitudes of 23, 25, 24, 23, and 22, respectively (McMahon 
et al. 2001). Magnitudes are put on a standart scale using 
observations of Landolt standard star fields taken on the 
same night. The accuracy of the preliminary photometric 
calibration is ±0.1 mag. 
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Figure 3. Colour-apparent magnitude diagram for the original 
sample. 



5.2 The overlapping sources, de-reddening of the 
magnitudes, bright stars, and extra-galactic 
objects 

The magnitudes are provided from the Cambridge As- 
tronomical Survey Unit (CASU). There are totally 14439 
sources in 24 sub-fields in the field SA 114. It turned out 
that 2428 of these sources are overlapped, i.e. their angular 
distances are less than \" to any other source. We omitted 
them, thus the sample reduced to 12011. The E(B - V) 
colour- excesses for the sample sources are evaluated by the 
procedure of Schegel, Finkbeiner, & Davis, (1998) and cor- 
rected by the following equations of Beers et al. (2002): 

E(B — V) = E(B - V)s, for (B - V)s < 0.10 (11) 
E(B -V) = 0.10 + 0.65[E(B - V) s - 0.10], 

for E(B-V)s >0.10 (12) 

where E(B — V) s is the colour-excess evaluated via the pro- 
cedure of Schlegel et al. The E(B — V) colour-excess con- 
tours for the field is given in Fig. 1 as a function of Galactic 
latitude and longitude. Then, the total absorption Ay is 
evaluated by means of the well known equation, 



Rv 



3.1 



(13) 



E(B - V) 

For Vega bands we used the R\/Rv data of Cox (2000) 
for the interpolation, where A = 358lA, 4846A, 6240A, 
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Table 2. The relation between the total absorptions for Vega 
bands and for V band of U BV photometry. 



Filter A e// (A) A e// (A) Ri/R v 



u 


3581 


638 


1.575 


24.3 


9 , 


4846 


1285 


1.164 


25.2 


r 


6240 


1347 


0.867 


24.5 


i 


7743 


1519 


0.648 


23.7 


z 


8763 


950 


0.512 


22.1 



7743A, and 8763A (Table 2), and derived R\ by their combi- 
nation with Ay. Finally, the de-reddened u , g , r , io, and 
z magnitudes were obtained from the original magnitudes 
and the corresponding R\. 

The histogram for the de-reddened apparent magnitude 
g (Fig. 2) and the colour- apparent magnitude diagram (Fig. 
3) show that there is a large number of saturated sources 
in our sample. Hence, we excluded sources brighter than 
g = 17. However, the two-colour diagrams (u — g )o - 
(9 —r )o and (g — r )o - (r — i )o in Fig. 4 indicate that there 
are also some extra-galactic objects, where most of them lie 
towards the blue as claimed by Chen et al. (2001). It seems 
that the star/extragalactic object separation based on the 
'stellarity parameter' as returned from the SExtractor rou- 
tines (Bertin & Arnouts, 1996) couldn't be sufficient. This 
parameter has a value between (highly extended), and 1 



(point source). The separation work very well to classify a 
point source with a value greater than 0.8. We adopted the 
simulations of Fan (1999), additional to the work cited above 
to remove the extragalactic objects in our field. Thus we re- 
jected the sources with (u — g )o < —0.10, and those which 
lie outside of the band concentrated by most of the sources. 
After the last process the number of sources in the sample, 
stars, reduced to 6418. The two-colour diagrams (u — g )o 
- (g — r )o and (g — r )o - (r —i)o for the final sample are 
given in Fig. 5. A few dozen of stars with (u — g )o ~ 0.10 
and (g — r )o ~ 0.20 are probably stars of spectral type A. 

5.3 Absolute magnitudes, distances, population 
types and density functions 

In the Sloan photometry, the blue stars in the range 15 < 
g < 18 are dominated by thick disk stars with a turnoff 
(9 ~ r ) ~ 0.33, and for g > 18, the Galactic halo, has 
a turnoff colour (g — r ) ~ 0.2, becomes significant. Red 
stars, (g — r ) > 1.3, are dominated by thin disk stars for 
all apparent magnitudes (Chen et al. 2001). In our case, 
the apparent magnitude which separates the thick disk and 
halo stars seems to be a bit fainter relative to the Sloan 
photometry, i.e. g ~ 19 (Fig. 3). Thus, stars bluer than 
(9 — r ) = 1.1 and brighter than g — 19 are separated 
to the thick disk population and the colour-magnitude dia- 
gram of 47 Tuc (Hesser et al. 1987) is used for their abso- 
lute magnitude determination, whereas those with the same 
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Figure 7. Spatial distribution for stars with absolute magnitude 6 < M(g ) < 7 as a function of apparent magnitude, (a) (17.0-17.5], 
(b) (17.5-18.0], (c)(18.0-18.5], (d)(18.5-19.0], (e)(19.0-19.5], (f)(19.5-20.0], (g) (20.0-20.5], (h)(20.5-21.0], and (i) (21-22] 
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Figure 6. (g — r )o colour histogram as a function of apparent 
magnitude, for the star sample, (a) for 17 < g Q < 22, (b) for 
17 < #0 — 19, and (c) for 19 < g' Q < 22. The vertical downward 
arrow shows the limit value (g — r )o = 1-1 mag which separates 
thin disk and thick disk - halo couple. 



colour but fainter than g = 19 are assumed to be as halo 
stars and their absolute magnitudes are determined via the 
colour- magnitude diagram of M92 (Stetson & Harris 1988). 
From the other hand, stars redder than (g — r ) = 1.1 are 
adopted as thin disk stars and their absolute magnitudes 
are evaluated by means of the colour-magnitude diagram of 
Lang (1992) for Pop I stars (Fig. 6). The distance to a star 
relative to the Sun is carried out by the following formula: 



M{g )) = 51ogr - 5 



(14) 



The vertical distance to the galactic plane (z) of a star could 
be evaluated by its distance r and its Galactic latitude (b) 
which could be provided by its right ascension and declina- 
tion. The precise separation of stars into different popula- 
tions have been carried out by their spatial distribution as 
a function of their absolute and apparent magnitudes. Fig. 
7 gives the spatial distribution of stars with 6 < M(g ) < 7, 
as an example, and Table 3 gives the full set of absolute and 
apparent magnitude intervals, and the efficiency regions of 
the populations. Halo stars dominate the absolutely bright 
intervals, thick disk stars indicate the intermediate and the 
thin disk stars the faint ones, as expected (Fig. 8). 

Number of stars as a function of distance r relative to 
the Sun, and the corresponding mean distance z* from the 
Galactic plane, for different absolute magnitude intervals 
for three populations are given in Tables 4-6. The logarith- 
mic space density functions, D* = logD + 10, evaluated by 
means of these data are not given in the tables for avoidance 
space consuming but they are presented in Fig. 9-11, where 
D = N/AVi,2, AVi, 2 = (7r/180) 2 (D/3)(r| - rf), □: size of 
the field (4.239 square-degree), n and T2'- the limiting dis- 
tance of the volume AVi,2, and N: number of stars (per unit 
absolute magnitude). The horizontal thick lines, in Tables 
4-6, corresponding to the limiting distance of completeness 
(zi) are evaluated by the following equations: 
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Table 4. Number of stars as a function of distance r relative to the Sun, and the corresponding mean distance from the Galactic 
plane, for different absolute magnitude intervals for thin disk (distances in kpc). Horizontal thick lines correspond the limiting distance 
of completeness. 



M(g') 


(5-6] 


(6-7] 


(7-* 


i] 


(* 


5-9] 


(9-10] 




(10-11] 


(11-12] 


(12-13] 


ri - r 2 


z* N 


z* N 


z* 


N 


z* 


N 


z* 


N 


z* 


N 


z* 


N 


z* N 


0.10-0.20 






















0.13 


30 


0.12 20 


0.20-0.30 














0.22 


4 


0.19 


42 


0.19 


115 


0.19 39 


0.30-0.40 














0.27 


21 


0.26 


83 


0.26 


129 


0.26 33 


0.40-0.60 










0.43 


59 


0.37 


66 


0.37 


154 


0.36 


182 


0.34 28 


0.60-0.80 






0.56 


14 


0.52 


65 


0.52 


58 


0.52 


158 


0.50 


71 




0.80-1.00 






0.67 


43 


0.68 


58 


0.68 


72 


0.67 


87 


0.66 


20 




1.00-1.25 




0.87 18 


0.84 


87 


0.82 


65 


0.85 


72 


0.83 


77 








1.25-1.50 




1.03 67 


1.02 


76 


1.03 


55 


1.03 


46 


1.01 


30 








1.50-1.75 
1.75-2.00 


1.26 8 
1.42 39 


1.21 60 
1.39 46 


1.21 
1.40 


59 
37 


1.20 
1.40 


66 
16 


1.20 
1.40 


37 
16 


1.24 


7 








2.00-2.50 


1.58 75 


1.71 41 


1.61 


57 


1.53 


9 


1.63 


9 












2.50-3.00 
3.00-3.50 


2.08 26 
2.50 11 


2.08 21 
2.24 2 










1.99 


4 












Total 


159 


255 




373 




393 


405 




638 




547 


120 



Table 3. Dominant regions for three populations; thin disk, thick 
disk, and halo, as a function of absolute and apparent magnitudes. 
The symbol z is the distance to the Galactic plane. 



M(g') 


9o 


Thin disk 


Thick disk 


Halo 


(12-13] 


(17-22] 


z < 0.50 






(11-12] 


(17-22] 


z < 0.80 






(10-11] 


(17-22] 


z < 1.30 






(9-10] 


(17-22] 


z < 1.50 


z > 1.50 




(8-9] 


(17-18] 


z < 0.60 


z > 0.60 






(18-19] 


z < 0.85 


z > 0.85 






(19-20] 


z < 1.25 


z > 1.25 






(20-22] 


z < 1.60 


z > 1.60 




(7-8] 


(17-18] 


z < 1.00 


z > 1.00 






(18-19] 


z < 1.30 


z > 1.30 






(19.0-19.5] 


z < 1.70 


z > 1.70 






(19.5-20.0] 




z < 2.55 


z > 2.55 




(20-22] 




z < 3.80 


z > 3.80 


(6-7] 


(17.0-17.5] 


z < 1.12 


z > 1.12 






(17.5-18.0] 


z < 1.40 


z > 1.40 






(18.0-18.5] 


z < 1.50 


z > 1.50 






(18.5-19.0] 


z < 1.90 


z > 1.90 






(19.0-19.5] 


z < 2.25 


z > 2.25 






(19.5-20.0] 




z < 4.60 






(20.0-20.5] 




z < 4.70 


z > 4.70 




(20.5-21.0] 




z < 5.60 


z > 5.60 




(21-22] 




z < 6.40 


z > 6.40 


(5-6] 


(17.0-18.0] 


z < 1.68 


z > 1.68 






(18.0-18.5] 


z < 2.18 


z > 2.18 






(18.5-19.0] 


z < 2.52 


z > 2.52 






(19.0-19.5] 




z < 4.50 


z > 4.50 




(19.5-20.0] 




z < 4.50 


z > 4.50 




(20-22] 






z > 4.50 


(4-5] 


(17-18] 




z < 2.20 


z > 2.20 




(18-19] 




z < 3.60 


z > 3.60 




(19-22] 






z > 5.00 



Table 5. Number of stars for different absolute magnitude inter- 
vals for thick disk (symbols as in Table 4). 



M(g) -> 


(4-5] 


(5-6] 


( 


6-7] 


(7-8] 


(8-9] 


r> — ri 


z* N 


z* N 


z* 


N 


z* 


N 


z* N 


0.5-1.0 














0.64 30 


1.0-1.5 


1.05 17 




1.13 


1 


1.04 


10 


0.95 59 


1.5-2.0 


1.29 51 




1.30 


95 


1.38 


40 


1.37 79 


2.0-2.5 


1.66 69 


1.78 80 


1.66 


136 


1.74 


65 


1.69 84 


2.5-3.0 


2.05 44 


2.05 116 


2.06 


124 


2.05 


99 


2.02 44 


3.0-3.5 


2.42 20 


2.42 115 


2.43 


129 


2.42 


67 


2.42 23 


3.5-4.0 


2.77 20 


2.81 107 


2.80 


104 


2.82 


53 


2.78 10 


4.0-4.5 


3.14 18 


3.18 63 


3.17 


72 


3.16 


47 


3.09 7 


4.5-5.0 


3.48 3 


3.57 66 


3.54 


69 


3.50 


30 


3.57 4 


5.0-5.5 




3.91 51 


3.90 


54 


3.75 


5 




5.5-6.0 




4.30 38 


4.26 


53 








6.0-6.5 




4.66 34 


4.62 


44 








6.5-7.0 




5.06 11 


5.03 


38 








7.0-8.0 




5.27 4 


5.38 


58 








8.0-9.0 






5.83 


24 








Total 


242 


685 


1001 




416 


340 




4 6 8 10 12 14 

M(g') 



(g -M(g)) = 5\ogri-5 (15) 
zi — risinb (16) 

where g is the limiting apparent magnitude (17 and 22, 
for the bright and faint stars, respectively), ri the limiting 
distance of completeness relative to the Sun, and M(g ) the 
appropriate absolute magnitude Mi or M2 for the absolute 
magnitude interval Mi — M2 considered. 



Figure 8. Absolute magnitude ranges dominated by different 
populations. Symbols: (+) halo, (□) thick disk, (•) thin disk, and 
thick line: the distribution of stars for all absolute magnitudes 

6 GALACTIC MODEL PARAMETERS 

6.1 Density laws for different populations and 

different absolute magnitudes: new approach 
for the model parameter estimation 

In the literature different density laws, such as secans hyper- 
hnlir.ns (sp.rh). snna.reH ser.a.ns Vivn^rhnlirns ( sp.rh 2 ) or fiynn- 
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Table 7. Galactic model parameters for different absolute magnitude intervals for thin disk resulting from the comparison of observed 
logarithmic space densities with a (unique) density law (Fig. 9). The columns give: M(g ) absolute magnitude interval, density law, 
logarithmic local space density n*, scale height for sech or sech 2 density law zo, scale height for exponential density law H, x 2 > standart 
deviation s, and local space densities for Hipparcos 0. 



M(g ) 


Density Law 


n* 


Zo/H (pc) 


2 /in-10\ 

x -(io ) 


s 





(12-13] 


exp 


^ .^-1-0 07 
8 - 40 ±0.07 


101 _7 


4362816 


± 0.14 


8.05 




sech 


8 14+ - 07 
° -0.07 


98/162 + J 7 , 


3767270 


0.14 






sech 2 


o n o + 0.06 


166/275 + ^2 


2691626 


0.12 




(11-12] 


exp 


o Qo + 0.02 


103 + ^ 


1940921 


0.10 


7.92 




sech 


8 69+ - 04 
a,D -0.04 


102/169+^ 


1690410 


0.09 






sech 2 


8 - 55 -0.03 


188/312+^ 


755904 


0.06 




(10-11] 


exp 


o qo + 0.04 
8 - 38 -0.04 


168 + ? 


628271 


0.10 


7.78 




sech 


o i -,+0.03 
8 - n -0.06 


168/278+^ 6 


857229 


0.10 






sech 2 


7 99+ - 04 
'• yy -0.04 


300/497+^ 


994846 


0.06 




(9-10] 


exp 


7 60+ - 05 
U -0.05 




274519 


0.10 


7.63 




sech 


7 3 4 +0.05 
' ,0 -0.05 


256/424+24 


295956 


0.06 






sech 2 


7 9q + 0.06 


460/762+^° 


370010 


0.06 




(8-9] 


exp 


7 47+0-07 
'•^ -0.06 


292+22 


247010 


0.13 


7.52 


(7-8] 


exp 


7 50+ - 02 
'• bU -0.02 


312+^ 


8207 


0.02 


7.48 


(6-7] 


exp 


7 43+0.02 
-0.02 


326+^ 


6295 


0.03 


7.47 


(5-6] 


exp 


7 43+0.01 

' -o.oi 


334+2 


162 


0.01 


7.47 



Table 6. Number of stars for different absolute magnitude inter- 
vals for halo (symbols as in Table 4). 



M(g) - 


(4-5] 




(5-6] 




(6-7] 




(7-J 


I] 


r\ — r 2 


z* 


N 


z* 


N 


z* 


N 


z* 


N 


2-4 


2.58 


19 










2.73 


14 


4-6 


3.81 


34 










4.17 


42 


6-8 


5.05 


25 


5.37 


99 


5.38 


76 


5.04 


48 


8-10 


6.57 


10 


6.71 


113 


6.65 


117 


6.72 


10 


10-15 


9.10 


27 


9.13 


178 


8.73 


129 






15-20 


12.97 


17 


12.54 


57 










20-25 


17.00 


3 


16.02 


2 










25-30 


19.44 


4 














30-35 


22.11 


2 














Total 




141 




449 




322 




114 



nentials, were used for parameterization of thin disk data, 
whereas the exponential density law was sufficient for thick 
disk data. In the present study, we compared the observed 
logarithmic space densities for thin disk with all density laws 
cited above. It turned out that sech 2 law fits better for three 
faint absolute magnitude intervals, i.e. 12 < M(g ) < 13, 
11 < M(g') < 12, and 10 < M(g) < 11, whereas exponen- 
tial law is favourable for brighter absolute magnitudes (Ta- 
ble 7 and Fig. 9). The argument used for this conclusion is 
the difference between the local space density resulting from 
the comparison of the observed space density functions with 
the density laws and the Hipparcos one (Jahreiss & Wielen 
1997). Table 7 gives also the corresponding scale- heights. It 
is interesting, the scale-height increase monotonously when 
one goes from the faint magnitudes towards the bright ones, 
however there is a discontunity at the transition region of 
two density laws. The range and the mean of the scale- height 
for thin disk are 264-497 pc, and 327 pc, respectively. Al- 
though 497 pc seems an extreme value, it is close to the 
upper limit cited by Robin & Creze (1986). 

For the thick disk, the observed logarithmic space den- 



Table 8. Galactic model parameters for thick disk. 712/ 'n\ indi- 
cates the local space density for thick disk relative to thin disk. 
Other symbols are same as in Table 7. 



M(g ) 


n* 


z /H (pc) 


x 2 -(io- 10 ) 


s 


n 2 /ni(%) 


(8-9] 
(7-8] 
(6-7] 
(5-6] 


6 19 + - 01 
D - -0.01 


970 + 27 
y -29 
806 + g 

895 + 2 8 3 

876+.f 5 


2002 
18505 
5879 
5503 


± 0.02 
0.07 
0.06 
0.08 


5.25 
6.46 
9.77 
9.55 



Table 9. Galactic model parameters for halo, hi and 723/721 give 
the axial ratio and the local space density for halo relative to thin 
disk, respectively. Other symbols are as in Table 7. 



M(g) 


n* 


K 


x 2 -(io- 10 ) 


s 


n 3 /ni(%) 


(7-8] 
(6-7] 


3.9812:1! 

4. a r£l» 


n 7R + 0.22 

°- 78 -0.18 

71+ - 10 
U - ' -0.08 


1959 
778 


± 0.30 
0.13 


0.02 
0.07 


(5-6] 
(4-5] 


4 53+ - 06 
4 79+O.O7 
4> ' y -0.06 


57+ - 03 
26+ - 01 


271 
189 


0.14 
0.16 


0.13 
0.31 



sity functions are compared with the exponential density 
law for the absolute magnitude intervals 8 < M(g ) < 9, 
7 < M(g') < 8, 6 < M(g') < 7, and 5 < M(g) < 6 (Table 
8, Fig. 10). The scale- height for different absolute magnitude 
intervals varies from 806 pc to 970 pc and are in agreement 
with the recent values in the cited literature (see the ref- 
erences in Table 1). The local space density relative to the 
one of thin disk, for the corresponding absolute magnitude 
intervals in Table 7, increases from the faint absolute mag- 
nitudes to the bright ones in a range (5.25-9.77)%, again 
in agreeable with the literature. The observed logarithmic 
space density functions for halo have been compared with 
the de Vaucouleurs density law for the absolute magnitude 



Galactic model parameters 9 




Figure 9. Comparison of the observed space density function with the density laws for different absolute magnitude intervals for the thin 
disk, (a) (5-6], (b) (6-7], (c) (7-8], (d) (8-9], (e) (9-10], (f) (10-11], (g) (11-12], and (h) (12-13]. Continuous curve: represents exponential 
law, whereas dashed curve, sech law and dashed-dot curve, sech 2 law. 



Table 10. Luminosity function resulting from the combination 
of local space densities for three populations, i.e. thin and thick 
disks, and halo, taken from Tables 7-9. The luminosity function 
of Hipparcos is given in the last column. 



M(g') 


<p*(M) 


s 





(12-13] 


8.08 


±0.14 


8.05 


(11-12] 


8.55 


0.10 


7.92 


(10-11] 


7.99 


0.10 


7.78 


(9-10] 


7.60 


0.10 


7.63 


(8-9] 


7.49 


0.13 


7.52 


(7-8] 


7.53 


0.02 


7.48 


(6-7] 


7.47 


0.03 


7.47 


(5-6] 


7.47 


0.01 


7.47 



intervals 7 < M(g ) < 8, 6 < M(g ) < 7, 5 < M(g ) < 6, 
and 4 < M(g) < 5 (Table 9 and Fig. 11). The axial ratio 
decreases from relative absolute magnitudes to the bright 
ones, whereas the local space density relative to the thin 
disk for the corresponding absolute magnitude intervals in- 
creases in the same order. The parameters cited here are 
in the range given in the literature, except the ones for the 
interval 4 < M{g ) < 5. The parameters derived for three 
populations have been tested by the luminosity function, 
given in Table 10 and Fig. 12, where ip*(M) is the total of 
the local space densities for three populations. There is a 
good agreement between our luminosity function and that 
of Hipparcos (Jahreiss & Wielen 1997) with the exception 
of absolute magnitude interval 11 < M(g ) < 12. Also, the 
corresponding standard deviations for all intervals are small 
(Table 10). 

We used the procedure of Phleps et al. (2000) for the 



1 1 1 1 1 1 1 1 1 1 

(a) 


i i i i i i i i i i 

(b) 


i i i i i i i i i i 


— — 

i i i i i i i i i i 



0123450123456 
z (kpc) z (kpc) 



I I I I 

(c) 


(d) 


[ I 





1 2 3 4 0.5 1 1.5 2 

z (kpc) z (kpc) 



Figure 10. Comparison of the observed space density function 
with the exponential density law for different absolute magnitude 
intervals for thick disk, (a) (5-6], (b) (6-7], (c) (7-8], and (d) (8-9]. 



error estimation in Tables 7-9 (above) and Tables 12 and 14 
(in the following sections), i.e. by changing the values of the 
parameters until % 2 increases or decreases by 1. 
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10 

z (kpc) 
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II 
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_ • _ 

I I I I I I I I I 


• 
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8 
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Figure 11. Comparison of the observed space density function 
with the de Vaucouleurs density law for different absolute mag- 
nitude intervals for halo, (a) (4-5], (b) (5-6], (c) (6-7], and (d) 
(7-8]. 




Figure 13. Comparison of the observed space density function 
combined for thin and thick disks and halo with the combined 
density law, for stars with 5 < M(g ) < 10. 





9 

M(g') 



10 



11 



12 13 



4 5 6 7 8 9 10 11 12 13 

M(g') 

Figure 12. Luminosity function obtained from combination of 
the local space densities for thin and thick disks and halo, result- 
ing from comparison of the observed space density function with 
the density laws, for different absolute magnitude intervals. (H) 
shows the Hipparcos values. 



6.2 Model parameter estimation by the procedure 
in situ 

We estimated the local space density and scale-height for 
thin disk and thick disk, and the local space density and 
axial ratio for halo simultaneously by the procedure in situ, 
i.e. by comparison of the combined observed space density 
functions with the combined density laws. We carried out 
this work for two sets of absolute magnitude intervals, 5 < 
M(g) < 10 and 5 < M(g) < 13. The second set covers 
the thin disk stars with 10 < M(g ) < 13 whose density 
functions behave differently than the density functions for 
stars with other absolute magnitudes. Number of stars as a 



Figure 14. Luminosity function resulting from the comparison of 
the combined observed space density function with the combined 
density law, for stars with 5 < M(g ) < 10. 



function of distance r relative to the Sun for eight absolute 
magnitude intervals are given in Table 11 and the density 
functions per unit absolute magnitude interval evaluated by 
these data are shown in Table 12 and Table 13 for the sets 
5 < M(g') < 10 and 5 < M(g) < 13, respectively. 



6.2.1 Model parameters by means of absolute magnitudes 
5 < M(g) < 10 

The observed space densities per absolute magnitude inter- 
val for three populations, i.e. thin and thick disks, and halo 
for stars with 5<M(g)<10, (Table 12) are compared with 
the combined density laws in situ (Fig. 13). The derived pa- 
rameters are given in Table 14. All these parameters are in 
agreement with the ones given in Table 1. Especially the 
relative local space density for thick disk, 8.32%, lies within 
the range given in two recent works (Chen et al. 2001 and 
Siegel et al. 2002). The resulting luminosity function (Fig. 
14) from the comparison of the model with these parameters 
and the combined observed density functions per absolute 
magnitude interval is also in agreement with the one of Hip- 
parcos (Jahreiss & Wielen 1997). However the error bars 
are longer than the ones in Fig. 12, particularly for the faint 
magnitudes. 
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Table 11. Number of stars as a function of distance r relative to the Sun, and the corresponding mean distance z* from the Galactic plane, 
for different absolute magnitude intervals (distances in kpc). Horizontal thick lines correspond the limiting distance of completeness. 



M(g') — 


(5-6] 




(( 


3-7] 


(7-1 


5] 


(8-9] 




(9-10] 


(10-11] 


(11-12] 


(12-13] 


r\ — r 2 


z* 


N 


z* 


N 


z* 


N 


z* 


N 


z* 


N 


z* 


N 


z* 


N 


z* N 


0.0-0.2 




















0.15 


1 


0.13 


30 


0.12 


20 


0.2-0.4 


















0.26 


25 


0.24 


125 


0.23 


244 


0.22 72 


0.4-0.7 














0.45 


77 


0.41 


96 


0.41 


236 


0.39 


230 


0.34 28 


0.7-1.0 










0.65 


57 


0.63 


135 


0.64 


100 


0.62 


163 


0.59 


41 




1.0-1.5 






1.00 


86 


0.93 


173 


0.93 


179 


0.92 


118 


0.89 


109 


0.75 


2 




1.5-2.0 


1.39 


47 


1.29 


201 


1.32 


136 


1.30 


149 


1.26 


53 


1.27 


4 








2.0-2.5 


1.69 


158 


1.68 


177 


1.68 


122 


1.67 


93 


1.63 


9 












2.5-3.0 
3.0-4.0 


2.06 
2.60 


142 
233 


2.06 
2.59 


145 
235 


2.05 
2.61 


99 
134 


2.02 
2.53 


44 
33 


1.99 


4 












4.0-5.0 


3.38 


129 


3.35 


141 


3.29 


77 


3.21 


10 
















5.0-7.5 
7.5-10.0 
10.0-12.5 
12.5-15.0 
15.0-20.0 


4.63 
6.53 
8.28 
10.21 
12.54 


209 
141 
100 
78 
57 


4.68 
6.40 
8.27 
9.90 


284 
181 
93 
36 


4.54 
6.42 


91 
14 


3.81 


1 
















20.0-25.0 


16.02 


2 




























Total 




1296 




1579 




903 




721 




405 




638 




547 


120 



Table 12. Logarithmic space density function, D* = log D + 10, 
per unit absolute magnitude interval for stars with 5 < M(g ) < 
10. < TV > is the weighted mean of number of stars for the ab- 
solute magnitude intervals in Table 11 under consideration, the 
other symbols are explained in the text (distances in kpc, volumes 
in pc 3 ). 



7*1 — 7*2 


AVi,2 


z* 


<N> 


D* 


0.4-0.7 


1.20 (5) 


0.408 


96 


6.90 


0.7-1.0 


2.83 (5) 


0.634 


118 


6.62 


1.0-1.5 


1.02 (6) 


0.929 


176 


6.24 


1.5-2.0 


1.99 (6) 


1.302 


162 


5.91 


2.0-2.5 


3.28 (6) 


1.678 


138 


5.62 


2.5-3.0 


4.90 (6) 


2.055 


129 


5.42 


3.0-4.0 


1.59 (7) 


2.600 


201 


5.10 


4.0-5.0 


2.63 (7) 


3.365 


135 


4.71 


5.0-7.5 


1.28 (8) 


4.656 


247 


4.29 


7.5-10.0 


2.49 (8) 


6.457 


161 


3.81 


10.0-12.5 


4.10 (8) 


8.279 


97 


3.37 


12.5-15.0 


6.12 (8) 


10.208 


78 


3.11 


15.0-20.0 


1.99 (9) 


12.536 


57 


2.46 



6.2.2 Model Parameters by means of absolute magnitudes 
5 < M(g) < 13 

We carried out the work cited in the previous paragraph 
for stars with a larger range of absolute magnitude, i.e. 
5 < M(g ) < 13. The observed density function is given 
in Table 13 and its comparison with the combined den- 
sity law is shown in Fig. 15. The derived parameters (Ta- 
ble 15), especially the local densities, are rather different 
than the ones cited in Sections (6.1) and (6.2.1). The reason 
for this discrepancy is that stars with absolute magnitudes 
10 < M(g ) < 13 have relatively larger local space densi- 
ties (Hipparcos; Jahreiss & Wielen 1997) and are closer to 
the Sun relative to stars brighter than M(g )=10, and they 
affect the combined density function considerably. Also the 
corresponding luminosity function is not in agreement with 
the one of Hipparcos, except one absolute magnitude inter- 
val, 12 < M(g ) < 13 (Fig. 16). 



Table 13. Logarithmic space density function, D* = log D + 10, 
per unit absolute magnitude interval for stars with 5 < M(g ) < 
13. < N > is the weighted mean of number of stars for the ab- 
solute magnitude intervals in Table 11 under consideration, the 
other symbols are explained in the text (distances in kpc, volumes 
in pc 3 ). 



7*1 — 7*2 


AVi, 2 




<N> 


D* 


0-0.2 


3.44 


(3) 


0.127 


25 


7.86 


0.2-0.4 


2.41 


(4) 


0.230 


147 


7.79 


0.4-0.7 


1.20 


(5) 


0.402 


187 


7.19 


0.7-1.0 


2.83 


(5) 


0.628 


133 


6.67 


1.0-1.5 


1.02 


(6) 


0.929 


176 


6.24 


1.5-2.0 


1.99 


(6) 


1.302 


162 


5.91 


2.0-2.5 


3.28 


(6) 


1.678 


138 


5.62 


2.5-3.0 


4.90 


(6) 


2.055 


129 


5.42 


3.0-4.0 


1.59 


(7) 


2.600 


201 


5.10 


4.0-5.0 


2.63 


(7) 


3.365 


135 


4.71 


5.0-7.5 


1.28 


(8) 


4.656 


247 


4.29 


7.5-10.0 


2.49 


(8) 


6.457 


161 


3.81 


10.0-12.5 


4.10 


(8) 


8.279 


97 


3.37 


12.5-15.0 


6.12 


(8) 


10.208 


78 


3.11 


15.0-20.0 


1.99 


(9) 


12.536 


57 


2.46 



Table 14. Galactic model parameters estimated by comparison of 
the logarithmic space density function for stars with 5 < M(g ) < 
10 (given in Table 12) and the combined density laws (Fig. 13). 
The symbols give: n*: logaritmic local space density, H: scale 
height, n/ny. local space density relative thin disk, and k= c/a: 
axial ratio for halo. 



Parameter 


Thin disk 


Thick disk 


Halo 


72* 


7 46+ - 02 

' -0.02 


6-38l° . 3 3 


4 30 +0 - 29 

^• OU -0.45 


if (pc) 


275±!j 


85lt£ 




n/m(%) 




8.32 


0.07 


K 






n ^7+0.33 

U -°' -0.23 



12 Karaali et al. 



Table 15. Galactic model parameters estimated by comparison of 
the logarithmic space density function for stars with 5 < M(g ) < 
13 (given in Table 13) and the combined density laws (Fig. 14). 
Symbols as in Table 14. 



Parameter 


Thin disk 


Thick disk 


Halo 


n* 


8.21 


6.22 


4.30 


H (pc) 


184 


1048 




n/m(%) 




1.00 


0.01 








0.67 





8t 

4 ^^^^^ 

3 — 

2 r * ■ ■ 

5 10 15 

z(kpc) 

Figure 15. Comparison of the observed space density function 
combined for thin and thick disks and halo with the combined 
density laws, for stars with 5 < M(g ) < 13. 



7 DISCUSSION 

We estimated the Galactic model parameters by means of 
the vertical density distribution for the field SA 114 (a = 
22> i 40 m 00 s ,<5 = 00°00 / 00 // ; I = 68°. 15, b = -48°.38; 4.239 
square- degree; epoch 2000) by means of two procedures, i.e. 
a new approach and the procedure in situ. The new ap- 
proach is based on the comparison of the observed space den- 
sity functions per absolute magnitude interval with a unique 
density law for each population individually, where the sep- 
aration of stars into different population types is carried out 
by a slight modification of the method of Karaali (1994), 
i.e. by their spatial distribution as a function of absolute 




6 

5 I i I i I i I i I i I i I i I i 

5 6 7 8 9 10 11 12 13 
M(g') 

Figure 16. Luminosity function resulting from the comparison of 
the combined observed space density function with the combined 
density law, for stars with 5 < M(g ) < 13. 



and apparent magnitude. This approach covers nine abso- 
lute magnitude intervals, i.e. 4 < M(g ) < 5, 5 < M(g ) < 6, 

6 < M(g') < 7, 7 < M(g') < 8, 8 < M(g') < 9, 
9 < M(g') < 10, 10 < M(g') < 11, 11 < M(g') < 12, and 
12 < M(g ) < 13. Whereas the procedure in situ compares 
the observed space density functions per absolute magni- 
tude interval with the combination of a set of density law 
representing all the populations. This procedure is carried 
out for two absolute magnitude intervals, 5 < M(g ) < 10 
and 5 < M(g ) < 13. However, we will not discuss the pa- 
rameters for 5<M(g)<13, since they are rather different 
than the parameters appeared in the literature up to now. 

7.1 Parameters determined by means of the new 
approach 

The new approach provides absolute magnitude dependent 
Galactic model parameters. The scale-height for thin disk 
increases monotonously from the faint magnitudes to bright 
ones, however there is a discontunity at the transition region 
of two density laws (exponential law for M(g ) < 10, and 
sech 2 law for 10 < M(g ) < 13). The range of this param- 
eter is 264-497 pc, and one can find all the values derived 
for different absolute magnitude intervals in the literature 
including the extreme one, i.e. 497 pc which is close to the 
one cited by Robin & Creze (1986). The thick disk is dom- 
inant in the intervals 5 < M(g) < 6, 6 < M(g') < 7, 

7 < M(g) < 8, and 8 < M(g) < 9. The scale height for 
thick disk lies within 806-970 pc and the range of the local 
space density relative to the thin disk for the corresponding 
absolute magnitude interval is (5.25-9.77)%, in agreement 
with the literature. The halo population is dominant in the 
intervals 4 < M(g) < 5 , 5 < M(g') < 6, 6 < M(g) < 7, 
and 7 < M(g ) < 8. The local space density and the axial 
ratio for the brightest interval are rather different than the 
ones in the other intervals and those cited up to now, this is 
probably due to less stars at large distances (see Table 6 and 
Table 9) . For the other absolute magnitude intervals the rel- 
ative local space density, ranging from 0.02% to %0.13, and 
the axial ratio with range 0.57-0.78 are also in agreement 
with the literature. The agreement of the luminosity func- 
tion, derived by the combination of the local space densities 
for three populations, with the Hipparcos one is confirma- 
tion for the Galactic model parameters. 

7.2 Parameters determined by means of the 
procedure in situ 

Let us reiterate the parameters derived by means of stars 
with 5 < M{g ) < 10. The scale-heights for thin and thick 
disks are 275 pc and 851 pc respectively, the axial ratio for 
halo is 0.67, and the local space densities for thick disk and 
halo, relative to the thin disk, are 8.32% and 0.07% respec- 
tively. These values are close to at least one, but not to all, of 
the values estimated for different absolute magnitude inter- 
vals by means of new approach (see Section 7.1 and Table 7). 
For example, the scale height 275 pc mentioned here is com- 
parable with the scale height for thin disk estimated for the 
interval 8 < M(g ) < 9, and all the parameters estimated 
via the procedure in situ are close to the corresponding ones 
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Table 16. The most appropriate Galactic model parameters 
(symbols as given in the previous tables). 





Thin 


disk 


Thick disk 




Halo 




M(g') 


H (pc) 


n l 


H (pc) 


n 2 /n 1 (%) 


K, 


n 3 /ni(%) 





(5-6] 


335 


7.4 


875 


9.5 


0.6 


0.15 


7.43 


(6-7] 


325 


7.4 


895 


9.8 


0.7 


0.05 


7.47 


(7-8] 


310 


7.5 


805 


6.5 


0.8 


0.02 


7.48 


(8-9] 


290 


7.5 


970 


5.2 






7.52 


(9-10] 


265 


7.6 










7.63 


(10-11] 


495 


8.0 










7.78 


(11-12] 


310 


8.6 










7.92 


(12-13] 


275 


8.1 










8.05 



estimated for the interval 8 < M{g ) < 9 by means of the 
new procedure. Also, these values are in agreement with the 
parameters cited in the literature. 

Although the errors are larger, the luminosity function 
resulting from the comparison of the observed space den- 
sity functions with the Galactic model with the parameters 
cited above agree with the luminosity function of Hipparcos 
(Jahreiss & Wielen 1997), confirming the parameters. 



the thin disk. This population is dominant in faint magni- 
tudes. For instance, since at late spectral types where stellar 
masses for the main-sequence stars are relatively small and 
as they occupy the space close to the Galactic plane, one 
would expect short scale heights for these absolute magni- 
tude intervals, whereas one expects large scale height for the 
ones corresponding to relatively bright absolute magnitude 
intervals. This is also the case in our work and is given in 
Table 7. When local space density for thin and thick disks, 
and for halo or for their both combination is considered, de- 
pendence of model parameters on absolute magnitude would 
be a result of the luminosity functions adopted up to now 
for individual populations. 

Some researchers restrict their work related with the 
Galactic model estimation with absolute magnitude. The 
recent work of Robin et al. (2003) who treated stars with 
M(V) < 8 can be given as an example. This is in confirmity 
with our work. Finally, we argue that the most appropriate 
Galactic model parameters are those estimated for individ- 
ual absolute magnitude intervals by comparison of the ob- 
served space density functions with the corresponding (and 
unique) density law. Our results obtained by this procedure 
are given in Table 16. 



7.3 How can we decide on the most appropriate 
Galactic model parameters? 

We estimated different sets of Galactic model parameters in 
two different ways, which both are in agreement with the 
ones appeared so far but differ from one another, when a 
specific one is considered. Now the question is: can we se- 
lect the most appropriate model parameters for our Galaxy? 
We emphasize that the model parameters are absolute - and 
hence mass - dependent. Also, absolutely faint stars do not 
contribute to thick disk and halo populations in the solar 
neighbourhood. Any procedure not regarding these two ar- 
guments will result in the estimation of model parameters 
with large range and very likely will differ with the pre- 
viously cited ones, depending on the absolute magnitude 
interval. 

The dependence of the model parameters on absolute 
magnitude can be explained as such; from the astrophysical 
point of view, the mass of a star and also its chemical compo- 
sition is two of the fundamental parameters during its forma- 
tion. Hence there is a good correlation between the mass of a 
main sequence star and its absolute magnitude, this in turn 
is also reflected in its spectral type. Therefore, dependence 
of a model parameters on absolute magnitude is expected. 
Particularly, if we consider the local space density of a thick 
disk population within the local volume, same number of 
stars with different magnitudes are not expected to exist. 
According to Fig. 8, the thick disk is mostly dominant by 
the stars in absolute magnitude interval of 6 < M(g ) < 7. 
If for different magnitude intervals, the same density law is 
used, then we can expect to get the highest relative local 
space density within the thick disk for the above given in- 
terval. Actually, the last column of Table 8 collaborates our 
expectation, where the highest relative local space density, 
which is 9.77% and the lowest relative local space density, 
which is 5.25%, correspond to absolute magnitude interval of 
6 < M(g ) < 7. The same holds for other model parameters 
when, for instance scale height is considered. Let us consider 
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